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Abstract. The (5-regime of QCD is characterised by light quarks in a small spatial box, but a 
large extent in (Euclidean) time. In this setting a specific variant of chiral perturbation theory 
— the 5-expansion — applies, based on a quantum mechanical treatment of the quasi one- 
dimensional system. In particular, for vanishing quark masses one obtains a residual pion mass 
M^, which has been computed to the third order in the <5-expansion. A comparison with nu- 
merical measurements of this residual mass allows for a new determination of some Low Energy 
Constants, which appear in the chiral Lagrangian. We first review the attempts to simulate 
2-flavour QCD directly in the 5-regime. This is very tedious, but results compatible with the 
predictions for have been obtained. Then we show that an extrapolation of pion masses 
measured in a larger volume towards the (5-regime leads to good agreement with the theoretical 
predictions. From those results, we also extract a value for the (controversial) sub-leading Low 
Energy Constant h. 
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1. QCD and Chiral Perturbation Theory 

The Lagrangian of QCD is scale invariant, but its quantisation singles out an intrinsic energy 
Aqcd 5 which sets the scale for the hadron spectrum. Our daily life is dominated by low energy 
and therefore by the lightest quark flavours, i.e. quarks with masses m q <C Aqcd- hi the limit of 
vanishing quark masses, their left- and right-handed spinor components l and ^ r) decouple. 
Thus the Lagrangian takes the structure 

£qcd = *lD^ l + $ R DV R + £ gauge , (1) 



where D is the Dirac operator. For Nf massless quark flavours, this Lagrangian has the global 
symmetry 

U(N f ) L U(N f ) R = SU(N f ) L <g> SU(N f ) R ® U(l) v <8> U(l) A . (2) 

chiral flavour symmetry baryon number conservation axial symmetry 

Here we split off the phases; the (vectorial) symmetry under simultaneous left- and right-handed 
phase rotation, J7(l)y, corresponds to the conservation of the baryon number. The remaining 
U(1)a symmetry — for opposite L and R phase rotations — is the axial symmetry. That is a 
symmetry of the classical theory, which breaks explicitly under quantisation, i.e. it is anomalous. 
We are interested in the remaining chiral flavour symmetry, which (in infinite volume) breaks 
spontaneously, 

SU(N f ) L ® SU(N f ) R — >• SU(N f ) L+R . (3) 

Chiral Perturbation Theory (xPT) deals with an effective Lagrangian in term of fields in the 
coset space of this spontaneous symmetry breaking, U{x) £ SU(Nf) p]. Thus it captures the 
lightest degrees of freedom, which dominate low energy physics, in this case given by Nj — 1 
Nambu-Goldstone bosons. The effective chiral Lagrangian C e g embraces all terms which are 
compatible with the symmetries. This concept also extends to the case where small quark 
masses are added, m q > 0. Then one deals with light pseudo Nambu-Goldstone bosons, which 
are identified with the light mesons; for Nf = 3 this includes the pions, the kaons and the 
77-meson. 

Here we consider the case Nf = 2, so we only deal with the quark flavours u and al. We 
assume them to be degenerate, i.e. to have both the mass m q . In this case the field U (x) € SU (2) 
describes the pion triplet. The terms in C e g are ordered according to an energy hierarchy, which 
depends on the number of derivatives and powers of m q . Some of the first terms are 

C cS = ^ 2 Tr[^C/t 5AtC /] + I Smq Tr[C/ + [/t] 

~h (Tr[d^d^U}) 2 - lz 2 (Tr[d^ d v U]f 

-03 + h)(^f + Urf + 1^ Tr[d,U%U] Tr[U + + (4) 

Each term comes with a coefficient, which is a free parameter of the effective theory. These 
coefficients are known as the Low Energy Constants (LECs). The leading LECs are the pion 
decay constant F n (which was measured as F n = 92.4 MeV) and the chiral condensate S 
(the order parameter for chiral symmetry breaking). The Zj are sub- leading LECs. In the 
chiral limit m q = 0, only F w occurs in the leading order. This determines an intrinsic scale 
A x = 4-7T F w — 1.2 GeV, which adopts the role of Aqcd in xPT. 

The LECs are of primary importance for low energy hadron physics. To some extent they can 
be fixed from phenomenology. On the theoretical side, they can only be determined from QCD 
as the underlying fundamental theory. Since this refers to low energy, it is a non-perturbative 
task, and therefore a challenge for lattice simulations. If one succeeds in their determination, 
we arrive at a rather complete, QCD-based formalism for low energy hadron physics. 

2. Pions in a finite volume: p-regime, e-regime and 5-regime 

For a system of pions in a finite volume, say with periodic boundary conditions and with a 
characteristic extent L, the low energy expansion can be formulated in terms of the dimensionless 
lowest non-zero momentum p^/K x ~ l/(2i ? 7r L), and of the correlation length £ = M~ l (M n is 



the pion mass). Both p^ and M v should be light compared to A x . Depending on the size and 
shape of the volume, one distinguishes various regimes, with different counting rules for these 
expansion parameters: 

• p-Regime: This is the standard setting with a large volume, L ^> £, and therefore small finite 
size effects [2J. In this case the counting is simply 0(1/ L) = 0(M n ). From the Lagrangian 
we can read off (on tree level) the Gell-Mann-Oakes-Renner relation 

Ml = m q . (5) 

• e- Regime: This regime refers to a small box in Euclidean space |3], say V = L 4 with 
L < £ (where £ is still the would-be inverse pion mass in a large volume, with the same 
quark mass; in Nature £ ~ 1.5 fm). In this regime, the xPT counting rules read 
0(1/ L) = 0(m q ) = 0(M%/A x ), unlike the p-regime counting. 

Experimentally the e-regime is not accessible, but QCD simulations in (or at least close to) 
this regime are feasible. They are of interest in particular because the LECs determined 
in this regime are the same that occur in large volume. Therefore we can extract physical 
information even from an unphysical regime. Since this can be achieved with a modest 
lattice size, this method is attractive from a practical perspective It has been intensively 
explored since 2003. It is difficult to simulate safely inside this regime, but certain 
properties, which are characteristic for the e-regime, have been observed. Ref. [5] provides 
a short overview. 

• 5-Regime: Here one deals with a small spatial volume, but a large extent T in Euclidean 
time [6], say 

L 3 xT , L<e«T . (6) 

This relation for L and T is depicted in Fig. [1] on the left (it is exactly opposite to the setting 
used in studies of QCD at finite temperature). The counting rules for the corresponding 
5 -expansion are 

A^I = °W ■ ^>A^ = °«- < 7 > 

A map of these three regimes in terms of the pion mass and the inverse extent in Euclidean 
time (the temperature) is shown in Fig. [I] on the right. 

This article addresses the <5-regime. It is far less known and explored than the p- and the e- 
regime, but it shares with the latter the exciting property that physical LECs can be extracted 
from an unphysical setting. A further motivation for studying QCD in a "<5-box" is that its 
shape allows (approximately) for a simplified analytical treatment in terms of Id field theory, 
i.e. quantum mechanics. In this case one considers a quantum rotator as described by the Id 
0(4) model, due to the local isomorphism between the orthogonal group 0(4) and the chiral 
symmetry group SU(2)l <g> SU(2)r. Closely related systems have applications in solid state 
physics, in particular regarding quantum anti-ferromagnets 0[8]. 

Spontaneous symmetry breaking does not occur in a finite volume. Therefore the piona3 - 
i.e. the pseudo Nambu-Goldstone bosons — cannot become massless in the chiral limit m q = 0, 
in contrast to the infinite volume. We may consider a fixed box of a J-shape and vary the quark 
mass: a large value of m q implies a large pion mass, so that we enter the p-regime and the 
Gell-Mann-Oakes-Renner relation ([5]), m q cx M%, is approximated. For small m q the pion mass 
turns into a plateau, which ends in the chiral limit at a residual pion mass M^. This behaviour 
is illustrated schematically in Fig. [2J 

1 One might argue if the term "pion" is adequate in the 5-regime. We find it acceptable and convenient, but 
readers who disagree may simply denote (see below) as the "mass gap" . 
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Figure 1. On the left: an illustration of a typical shape of a 5-box, i.e. an anisotropic finite 
volume where a pion gas can be treated by the (^-expansion. On the right: a schematic map 
of the applicability domains of three different expansion rules of xPT, namely the p- , the e- 
and the ^-regime. The dashed lines indicate regions where clearly one expansion holds; in the 
transition zones between these regions various expansions could work more or less. 
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Figure 2. A qualitative picture of the expected behaviour of the pion mass squared in a <5-box. 
For heavy quarks and pions we approximate the p-regime relation m q oc M%. For light quarks 
the pion mass attains a plateau, and finally (in the chiral limit m q = 0) the residual value M^. 



The value of M* can be computed with the ^-expansion. The spectrum of the 0(4) quantum 
rotator (a quantum mechanical particle on the sphere S* 3 ) is given by E% = l{£ + 2)/(20), so the 
mass gap amounts to = 3/(20). The challenge is now to compute the moment of inertia 
O. In his seminal paper on the 5-regime, H. Leutwyler gave its value to leading order (LO) as 
O F%L 3 . Thus the residual pion mass can be written as 



The shift A captures higher order corrections, which are suppressed in powers of 1/(F W L) 2 . 
They have been evaluated to next-to-leading order (NLO) in Ref. [8j, and recently even to 



next-to-next-to-leading order (NNLO) [9], which yields 



0.4516... 0.08843... 
A = ^ OTO + 



1 -0.1599... (ln(AiL) + 41n(A 2 L)) . (9) 



Aj are scale parameters for the sub-leading LECs. The latter are given at the scale of the 
physical pion mass as 

U = ln(Aj/MP hys ) 2 . (10) 

Even more recent papers addressed again the NNLO of the (^-expansion [10] . and the corrections 
due to finite 7n q [11] , In the following we will discuss numerical results for M^. We see that a 
confrontation with the analytical prediction in eqs. ([5]), Q could enable a new determination of 
a set of LECs from first principles of QCD. 



3. Attempts to simulate QCD in the <5-regime 

The straight way to measure are simulations directly in the 5-regime. Since the 5-box 
differs from the lattice shapes in usual simulations, this requires the special purpose generation 
of configurations. Moreover, precise chirality is vital in this regime, hence one is supposed to 
use a formulation of lattice quarks which preserves chiral symmetry. Such lattice fermions are 
known since the late 90ies, but their simulation is extremely tedious, in particular with dynamical 
quarks (i.e. keeping track of the fermion determinant in the generation of gauge configurations). 
We anticipate that so far there are no robust results of simulations clearly inside the 5-regime. 
In this section we summarise the efforts that have been carried out so far. 

At the Symposium LATTICE2005 D. Hierl presented a first attempt to simulate 2-flavour 
QCD in the 5-regime |12j . That study used a truncated version of a chiral lattice Dirac operator, 
so a first question is if the quality of approximate chirality was sufficient for that purpose. Since 
this Dirac operator is very complicated, that simulation was performed with a non-standard 
algorithm, which probes the fermion determinant with a stochastic estimator. The spatial 
volume was ~ (1.2 fm) 3 , and the results for M w at small quark masses agreed well with the LO 
of the ^-expansion, i.e. eq. ([8]) at A = (and with the phenomenological value of F n ). 

In 2007 the QCDSF Collaboration generated a new set of data, which have not been published. 
They were obtained with dynamical overlap quarks, which are exactly chiral. This simulation 
used the Hybrid Monte Carlo algorithm (i.e. the reliable standard algorithm). The lattice had 
a modest size of 8 3 x 16 sites, and the spatial box length was again L ~ 1.2 fm. At first sight 
the results seemed to look fine: for decreasing m q we saw a transition from a Gell-Mann-Oakes- 
Renner type behaviour towards a plateau. Its value agreed with the chirally extrapolated value 
of Ref. [12], and therefore also with the LO of eq. (|8|). 

Unfortunately, this is not the end of the story. If we proceed to the NLO correction, i.e. if 
we include the first term of A given in eq. (j9]) , the predicted value for decreases drastically 
in this small volume — from 782.5 MeV down to 321.7 MeV — and the agreement with the 
above data is gone. Considering this dramatic effect of the NLO correction one might worry 
that the <5-expansion could converge only very slowly in this small box, and such simulations are 
not instructive at all. However, adding also the NNLO correction alters the NLO result only a 
little — to (336.3 ± 7.6) MeV — so it is reasonable to assume the ^-expansion to be already well 
converged o 

Thus simulations in this small box can be useful in view of a confrontation with the analytical 
predictions for M^. In fact a second sequence of runs by the QCDSF Collaboration (performed 

2 To obtain this theoretically predicted value, we inserted the LECs as far as they are known. In particular the 
sub-leading LECs li = —0.4 ±0.6 and h = 4.1 ±0.1 are taken from Ref. [13] . along with their uncertainties, which 
imply the uncertainty in the NNLO value of M^- 
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Table 1. The results by QCDSF Collaboration of the year 2008, on a 8 3 x 16 lattice with 5 
values of the mass m q for the two degenerate dynamical overlap quark flavours. We display 
the measured values of the physical lattice spacing a and the pion mass, as well as the NNLO 
prediction for the residual pion mass M*. (The errors capture the uncertainty in a, in the pion 
mass in lattice units, and in the LECs l\ and 1%, cf. eq. (|14p .) 

in 2008) yielded data much closer to the NNLO prediction. They involved 5 quark masses; 
the results for the physical lattice spacing a and the pion mass are given in Table [TJ Since a 
varies for the different values of m q , the box length L = 8a and the prediction for vary as 
well. In Fig. Owe compare the numerically measured pion masses and the corresponding NNLO 
^-expansion results. They are compatible, in contrast to the earlier data, which were probably 
not well thermalised. In these new runs thermalisation is accomplished, but the ratio T/L = 2 
is still modest. We therefore proceeded to a 8 3 x 32 lattice, where our runs are ongoing. Still, 
the quality of agreement with the <5-expansion that we observed already on the 8 3 x 16 lattice 
is impressive. 
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Figure 3. An illustration of our results given in Table [1} we show the numerically measured 
pion mass and the theoretical residual mass M^, as predicted by the NNLO ((-expansion. (The 
latter actually refer to the chiral limit, but its value depends on L and therefore on the lattice 
spacing a, which varies for the simulations at different m q .) We observe good agreement between 
the measured and the predicted values. 



4. Residual pion mass by an extrapolation from the p-regime 



In this section we proceed to a different approach. It is based on simulation results in the p- 
regime (up to the transition zone), which are then extrapolated towards the J-regime. Details of 
this study are given in Ref . [2] . Also in this framework we consider it essential to use dynamical 



quarks, but we do not insist on exact chiral symmetry in the p-regime. Hence we used Wilson 
fermions, which is an established standard lattice fermion formulation, in a form which corrects 
0(a) lattice artifacts. Thus the simulation was much faster, and we could tackle much larger 
lattices than those mentioned in the Section 3; our data reported below were obtained on three 
lattice sizes: 24 3 x 48 , 32 3 x 64 and 40 3 x 64 . On the other hand, this lattice regularisation breaks 
the chiral symmetry explicitly, so that additive mass renormalisation sets in. Nevertheless we 
were able to attain very light pion masses. 

For the gauge part we used the standard plaquette lattice action. Our simulations were 
carried out at two values for the strong gauge coupling g s , respectively the parameter (3 = Q/g 2 - 
We determined the physical lattice spacing a from the measured nucleon mass, which revealed 
that we were dealing with fine lattices, 

P = 5.29 -> a ~ 0.075 fm , (3 = 5.4 -> a ~ 0.067 fm . (11) 

Thus the spatial size was in the range L ~ 1.6 fm...3.0 fm. Due to the additive mass 
renormalisation, we could not refer to the bare quark mass anymore. We measured the current 
quark mass by means of the PCAC relation, 



m, 



(d 4 Ai(0,x A )P(0)) 
' <P(0> 4 )P(0)) ' 



where P is the pseudoscalar density, and is the axial current. We observed practically no 
finite size effects on a and on m q , but we did see a striking L-dependence of M n , as expected. 
These quantities were found in the range 

m q = 3.60 MeV ... 231 MeV , M ff = 174 MeV ... 1.52 GeV , M^L = 2.7 ... 9.7 . (13) 

The latter confirms that our data range from the deep p-regime to the transition zone. A few 
missing data points have been completed with a (lengthy) formula for exponentially suppressed 
finite size effects ]Yb\ I16| within the p-regime, which holds up to 0(p 4 ). This formula involves 
the renormalised sub-leading LECs k, i = 1 . . . 4, see eqs. (j3D and (fTU|h Well established 
phenomenological values were composed in Ref. [J3] (li, I2 were anticipated in footnote [2]) , 

[ 1 = -0.4 ±0.6 , F 2 = 4.1±0.1 , [ 3 = 2.9 ±2.4 , /" 4 = 4.4±0.4 . (14) 

They were estimated from irir scattering data, and in particular the I4 value is based on the 
scalar pion form factor. The dark horse in this context is Z3: we replaced the above value by 
Z3 ~ 4.2, which we obtained in our study, see below. 

Our measured and interpolated data are given in Ref. [13] ■ We extrapolated them towards 
the (5-regime and extracted in particular a value for based on the relatively simple chiral 
extrapolation formula 

M n (L) 2 = M* :2 + Cim q [l + C 2 m q ln(C 3 m q )] , (15) 

which interpolates between the 0(p 4 ) correction formula (for large LM V ) and the chiral limit 
[13] . The C{ and are treated as free parameters to be fixed by the fit. Our data and the fits 
for (3 = 5.29 and for /3 = 5.4 are shown in Figs. H]and[5l respectively. These plots also illustrate 
the extrapolation result for in the chiral limit. 

Our main result is shown in Fig. [6) It compares the extrapolated values for with the 
predictions based on the 5-expansion to LO, NLO and NNLO, as a function of L. The latter two 
predictions are very close to each other for the volumes under consideration, so we can again 
assume the expansion to be well converged. 




Figure 4. Our chiral extrapolation referring to the lattice spacing a ~ 0.075 fm (corresponding 
to (3 = 5.29) with L ~ 1.8 fm (on the left) and L ~ 3.0 fm (on the right). The data points 
show the measured pion mass squared against the quark mass in the p-regime (the Sommer scale 
parameter ro = 0.467 fm is employed to convert them into dimensionless units). The curve is 
the fit according to eq. (|15p . It ends in the chiral limit, where we illustrate the extrapolated 
value for and its error. 




Figure 5. The same as Fig. [U but now for the finer lattices with a ~ 0.067 fm (corresponding 
to /3 = 5.4), such that L ~ 1.6 fm (on the left) and L ~ 2.1 fm (on the right). 

The extrapolation results for reach down to values even below the physical pion mass. 
The plot also shows the line where the product M^L decreases below 1; this can be roughly 
considered as the boundary of the 5-regime. Our extrapolations lead close to this boundary, and 
in one case into the <5-regime|f] 

In particular this plot shows that the extrapolated masses match this curve remarkably 
well. It would not have been obvious to predict this feature, because our data were obtained in 

3 In that case, the extrapolation should actually turn into a plateau at tiny m q , but this would hardly change 
the result. 
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Figure 6. A comparison of our extrapolated values for the residual pion mass with the 
predictions by the 5-expansion (LO, NLO, NNLO), as functions of the spatial box size L. We 
observe a remarkably good agreement. In our largest volume, the extrapolated is located 
below the physical pion mass, and in the region M*L < 1, which can be viewed as the domain 
of the (5-regime. 

a regime where the basis of the ^-expansion (frozen spatial degrees of freedom) does not hold. 
This suggests that the extrapolation formula (fl5|) applies well in a sizable domain. 

5. Evaluation of Low Energy Constants 

By fitting our pion mass data according to eq. (fl~5j) . we obtained results for the four free fitting 
parameters. In Section 4 we discussed the results for that we obtained in this way. Moreover, 
the results for Ci, C<i and C3 can be used to evaluate the notorious sub-leading LEC £3 [2]. 
This is how we obtained our value that we anticipated in Section 4, 

h = 32ir 2 F% ^ In C \ = = 4.2 ± 0.2 . (16) 

(The error given here emerges from the fits, the additional systematic error would be hard to 
estimate). In view of other results in the literature, this value is in the upper region. An overview 
has been presented in Ref. [T7], in particular in Table 11 and Figure 9. That overview estimates 
the world average as I3 = 3.3(7). 

In the previous considerations of Sections 4 and 5, we inserted the phenomenological value 
of the pion decay constant, F n = 92.4 MeV. Alternatively, we could also treat F n as a free 
parameter to be determined by the fits. In particular, matching this yields in the chiral 
limit [Ti] 

^numerical = 7g +14 Me y (17) 

m q =0 

which seems a bit low. However, effective field theory considerations suggest that the value 
of F^ in the chiral limit should indeed be below the physical value; Ref. [IT] estimates 
FJ m -o ~ 86 MeV. 



6. Conclusions 

The (5-regime refers to a system of pions in a finite box, typically of the shape L 3 x T with 
L <M~ l <C T. It can be treated by Chiral Perturbation Theory with suitable counting rules, 
the 5-expansion. In particular this predicts the residual pion mass in the chiral limit, M^, which 
has been computed recently to NNLO [9]. 

First attempts to simulate 2-flavour QCD in the 5-regime were confronted with technical 
difficulties, in particular thermalisation problems. Nevertheless we obtained good results for the 
pion mass with light quarks on a small lattice of size 8 3 x 16, with L ~ (0.82 . . . 0.91) fm. 

In another pilot study we measured pion masses in the p-regime (up to the transition region) 
and extrapolated them towards the (5-regime. This yields numerical results for M^, which 
agree remarkably well with the predictions by the ^-expansion. That expansion is based on 
assumptions, which do not hold in the regime where the data were obtained. Therefore it would 
not have been obvious to predict the observed agreement. This comparison can be viewed as a 
numerical experiment, which led to an interesting observation. 

Our pion mass fits from the p- towards the <5-regime also fix some further constants, which 
allow for a new determination of the (mysterious) sub- leading Low Energy Constant I3] we 
obtained I3 = 4.2(2). This is somewhat above the average of the phenomenological and numerical 
estimates in the literature. The constants h,---,li are relevant in the effective description of tttt 
scattering. Finally we also considered the option to treat F n as a free parameter. In the chiral 
limit the fits lead to a value somewhat below the phenomenological F n . 

Robust numerical results, measured manifestly inside the 5-regime, are still outstanding. 
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